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ABSTRACT 

We investigate the behavior and consequences of the reverse shock that terminates the supersonic expansion of the baryonic wind 
which is driven by neutrino heating off the surface of (non-magnetized) new-born neutron stars in supernova cores. To this end we 
perform long-time hydrodynamic simulations in spherical symmetry. In agreement with previous relativistic wind studies, we find 
that the neutrino-driven outflow accelerates to supersonic velocities and in case of a compact, ~ 1 .4 M Q (gravitational mass) neutron 
star with a radius of about 10 km, the wind reaches entropies of about 100 per nucleon. The wind, however, is strongly influenced 
by the environment of the supernova core. It is decelerated and shock-heated abruptly by a termination shock that forms when the 
supersonic outflow collides with the slower preceding supernova ejecta. The radial position of this reverse shock varies with time and 
depends on the strength of the neutrino wind and the explosion conditions in progenitor stars with different masses and structure. Its 
basic properties and behavior can be understood by simple analytic considerations. We demonstrate that the entropy of the matter 
going through the reverse shock can increase to a multiple of the asymptotic wind value. Seconds after the onset of the explosion it 
therefore can exceed 400 k s per nucleon in low-mass progenitors around 10M Q , where the supernova shock and the reverse shock 
propagate outward quickly. The temperature of the shocked wind has typically dropped to about or less than 10 9 K, and density and 
temperature in the shock-decelerated matter continue to decrease only very slowly. For more massive progenitors with bigger and 
denser metal cores, the explosion expands more slowly so that the termination shock stays at smaller radii and affects the wind at 
higher temperatures and densities. In this case the termination shock might play a non-negligible, strongly time- and progenitor- 
dependent role in discussing supernova nucleosynthesis. 

Key words, supernovae: general — neutrinos — nuclear reactions, nucleosynthesis, abundances — hydrodynamics 

Neutron stars are born as extremely hot and dense remnants driven wind", unavoidably accompanies the birth of a hot, 
at the center of exploding massive stars. Shortly after their for- neutrino-cooling neutron star, independent of the details of the 
mation they heat up to temperatures that can reach tens of MeV not finally understood supernova explosion mechanism (in case 
and become ev en higher than 50 MeV for soft nuclear equations of sufficiently strong magnetic fields, however, there my be mag- 
^ , of state ( e -g- | Burrows & Lattimetj Il986t iKeil & Jankal 119951; netically drive n outflow instead of a thermally-driven wind; see 
- - 1 iPons et all 1 19991) . Their gravitational binding energy is carried iMetzger et al.l l2006). The mass loss of the nascent neutron star 
away by neutrinos, which are abundantly produced at such con- begins after the supernova explosion has been launched and con- 
ditions. On a timescale of seconds these neutrinos diffuse out of tinues until the neutron star is essentially transparent to neutri- 
the interior of the star and escape from their mean surface of last nos. This flow of baryonic matter is a rapidly expanding and 
scattering, the neutrinosphere. By the associated energy and lep- cooling high-entropy environment, in which an a-rich freeze- 
ton number loss, the hot, still proton-rich and neutrino-opaque out can lead to the production of elements heavier than the iron 
proto-neutron star thus evolves to the final cold, neutron-rich and group. For sufficiently large neutron exce ss, the g-rich freeze- 
neutrino-transparent remnant during roughly the first minute of out can merge smoothly into an r-process ( Woo slev & Hoff man 
its life. ll992HMeveret al.lll992l) . 

Outside of the neutrinosphere the radiated neutrinos, which 

have typical energies of 10-20 MeV, travel through a layer with In fact < a number of parameters has been recognized to 

a very steep density gradient and decreasing temperature. It is determine the possibility of r-process nucleosynthesis in the 

unavoidable that the residual interactions of the high-energy neutrino-wind environment: The neutron-to-proton ratio in the 

neutrinos with the cooler stellar matter deposit energy in this wind < expressed in terms of the electron-to-baryon ratio or 

region. This energy transfer does not allow the "surface" lay- electron faction Y e ; the expansion timescale, t, which decides 

ers of the hot, neutrino-cooling neutron star to remain in hy- how fast me temperature and density of the outflowing mat- 

drostatic equilibrium, but leads to mass lo ss at a low rate in a ter dro P; and the wind entropy per nucleon, s as a measure of 

neutrino-driven outflow of baryonic matter dDuncan et alll98l the photon-to-baryon ratio of the environme nt flWitti etalj|l994 

IWooslev & Ba7ol3[T992h . This outflow, the so-called "neutrino- iQian & Woosley||l996fc|Hoffman et alj|l996|). The entropy of the 

wind is typically tens to more than 100 per nucleon, making 

Send offprint requests to: the wind environment a candidate for the so-called high-entropy 

Correspondence to: thj@mpa-garching.mpg.de r-process dMever et al.ll 9921: iMeveJ 19 94). In addition, the mass 
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loss rate decides whether the wind could be the major source 
of the observed galactic abundances of r-process material. These 
wind parameters depend on the neutron star properties, in partic- 
ular on the gravitational field of the neutron star and thus its mass 
and radius, and on the neutrino emission of the neutron star, i.e., 
the ti me-dependent luminositi es and spectra of the radiated neu- 
trinos dOian & Woos lev 1996). Since it is mainly the absorption 
of electron neutrinos, v e , and antineutrinos, v e , on free neutrons 
and protons, respectively, which heats the stellar gas and is re- 
sponsible for driving the mass loss and for setting the electron 
fraction in the ejected gas, the e mission proper ties of these neu- 
trinos are most relevant (but see Wanaiol l2006l) for a suggestion 
that additional energy input to the wind by vv-annihilation in 
case of largely anisotropic neutrino emission from the nascent 
neutron star could lead to a decisive increase of the wind en- 
tropy). 

Only in case the neutrino-driven outflow becomes super- 
sonic beyond a critical point, the sonic point, it truly deserves 
the name "wind". In such wind solutions the physical con- 
ditions at the neutrinosphere and behind the supernova shock 
are causally disconnected. The presence of the sonic point un- 
ambiguously determines the solution for a given value of the 
driving luminosity. Wind solutions possess the highest ("criti- 
cal") mass loss rate (and the lowest specific total energy of the 
ejected matter) for a given neutrino luminosity. Physical solu- 
tions with larger mass loss rates (and lower specific total en- 
ergy) do not exist. Lower mass loss rates ( higher specific total 
energies) correspond to "breeze solutions" dTakahashi & Jankal 
1 1 997t [Otsuki et al1l2000l) . In these, the outflow velocity reaches 
a maximum and then decreases again to asymptote to zero at 
infinity. The whole region between the proto-neutron star sur- 
face and the outer boundary of the considered outflow is there- 
fore in sonic contact. While wind solutions are characterized by 
a continuously rising velocity and decreasing temperature, the 
temperature of breezes level off to a constant value at large radii 
where the flow is dominated by internal instead of kinetic en- 
ergy. This limiting value of the temperature at large distances 
from the neutron star is an additional characteristic parameter of 
breeze solutions. 

Transsonic neutrino-driven winds in the context of super- 
nova explosions and nucleosynt hesis were investigated by means 
of hydrodynam ic simulation s (Wooslev & Baron 1992), ana - 
lytic discussion dOian & Wooslevlll996l:ICardall & Fuller! 19971) . 
and numerical solutions of the steady-s t ate wind equations 
dOtsuki et al.ll2000t iThompson etai]|200ll) . lOtsuki et al.1 d~2000) 
discussed the differ ence between winds and breezes, but like 
IWanaio et al.l d200ll) they concentrated on the subsonic solutions 
for their nucleosynthesis calculations, mainly because these al- 
lowed them to set a boundary value of the temperature at some 
large radius. This was understood to mimic the transition of the 
wind into a dense shell of ejecta behind the outgoing supernova 
shock, the presence of which hampered the free expansion of the 
wind. Such a behavior was found in calculations of supernova 
explosions by the Li vermore group, which were e mployed in the 
r-proc ess st udies of IWooslev & Hoffman ( 1992), IWooslev et al.l 
(1994j), and iHoffman et all d 1996b . and als o in hydrodynamic 
sim ulations of neutrino-d riven outflows by IWitti et aTl (l994) 
and|Takahashi et alJd 19941) . which were started from postbounce 
models provided by the Livermore group. The outflow trajec- 
tories in these simulations showed temperature and density de- 
clining asymptotically to nearly constant values, which were 
reached when the flow was gradually decelerated upon catching 
up with the slower, ea rlier ejecta behind the supe rnova shock. 
Sumivosh i et al.l d2000t) and lTerasawa et all d2002l) also referred 



to this behavior for using an artificially imposed constant pres- 
sure at the outer boundary in their Lagrangian hydrodynamic 
simulations of neutrino-driven mass ejection. The external pres- 
sure produced outflow deceleration similar to that found in the 
previous supernova models. 

Applying modern, high-resolution shock-capturing schemes 
and a better numerical resolution to l ong-time hydrodynamic 
simul ations of sup e rnova explosions, iJanka & Miiller ( 1995b 
and iBurrows et al. d 19951) (see also the more recent mod- 
els of Bur as et al. 2006b used for nucleosynthesis studies by 
iPruet et ail 120051) discovered the formation of a wind termi- 
nation shock caused by the collision of a transsonic neutrino- 
driven wind with the dense, slower ejecta shell behind the su- 
pernova shock. So far, however, this reverse shock, which leads 
to an abrupt deceleration and shock heating of the fast wind, 
has not received much attention. Subtle, potentially significant 
effects in the r-process nucleosynthesis that may depend in in- 
teresting ways on the location of and streng th of the reverse 
shock were found by Thompson et al.l d2001l) . Although these 
authors mentioned a rather modest reheating of the wind mate- 
rial by the reverse shock passage (that causes a increase of the 
specific entropy of 10 per nucleon), they obtained a consid- 
erably enhanced production of third-peak r-process nuclei due 
to a slower postshock expansion and a significantly higher tem- 
perature (0.05 MeV instead of 0.01 MeV for unsh ocked winds) 
at the time the r-process freeze-out happens. Also IWanaio et al.l 
(2002), alluding to the possibility of a wind termination shock, 
introduced a freeze-out value TV as the final temperature of the 
wind, i.e., they limited the temperature (and density) decrease 
in the supersonic wind by a chosen lower value. The choice of 
this temperature was , naturally, to some degree ad hoc, although 
Wanaio et al. (2002) justified it by nucleosynthesis considera- 
tions. A systematic and detailed exploration of the formation of 
the wind termination shock, of its hydrodynamic al effects on the 
wind properties, and of its nucleosynthetic consequences, how- 
ever, is still lacking. 

In this paper we will study the time-dependent evolution of 
the wind termination shock in different progenitors with spheri- 
cally symmetric (ID) models. For this purpose we perform sim- 
ulations of neutrino-driven explosions, employing the approx- 
imations to the full supe rnova physics used in previous works 
dScheck et al.l2 004, 2006). The neutron star in our simulations is 
replaced by a contracting inner boundary at which neutrino lumi- 
nosities are imposed such that supernova explosions with a typ- 
ical explosion energy of 1-2 x 10 51 erg = l-2bethe (B) are trig- 
gered by neutrino heating. The subsequent explosion and evo- 
lution of the relic neutron star is followed until 10 s after core 
bounce. Varying the neutron star contraction, which depends on 
the incompletely known high-density equation of state, and the 
time-dependent neutrino emission from the forming neutron star, 
we will also investigate the sensitivity of the reverse shock ef- 
fects on the neutrino-wind properties. Our results suggest that 
wind termination shocks are a robust, long-lasting feature in the 
supernova core just like the outgoing supernova shock and the 
neutrino-driven wind are. Of course, since a final understand- 
ing of the explosion me chanism of core-collapse supernovae 
is still missing (see, e.g. iBuras et al.ll2003l l2006bllal and refer- 
ences therein) and because we excise the neutron star at the grid 
center instead of simulating its neutrino-cooling evolution and, 
moreover, make severe approximations to the imp ortant neu- 
trino transport (see Sect. 2.2 and Appendix D in IScheck et al.l 
2006), our calculations will not be able to yield final answers. 
Nevertheless, our results are suitable for discussing fundamen- 
tal properties of the wind termination shock and for developing 
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a basic understanding of the role of this so far not well studied 
aspect of supernova explosions. 

Our paper is structured in the following way. In Sect. Q] we 
will briefly describe the numerical approach taken in this work. 
In Sect. [2] we will present some comparisons we made with pre- 
vious work, in particular testing our approximative treatment of 
relativistic effects. Section [3] contains our results for spherically 
symmetric simulations, providing a reference case and then in- 
vestigating varied conditions (neutron star radius, neutrino lu- 
minosities) at the lower grid boundary and different progenitor 
stars. In addition, we will present an analytic discussion which 
allows one to basically understand the behavior and the proper- 
ties of the wind termination shock. Section [5] will finish with a 
summary and conclusions. 

1. Numerical aspects 

1.1. Hydrodynamics and neutrino treatment 

The simulations of this paper were carried out with the 
ne utrino-hydrod y namic s code and the microphysics described 
by IScheck et al.l (I2006I) . The hydrodynamics module is a ver- 
sion of the Prometheus code which is based on a direct Eulerian 
implementation of the Piec ewise Parabolic Method (PPM) of 
IColella & Woodward! (I1984 . It is a high-resolution shock cap- 
turing scheme and performs a conservative, explicit integration 
of the Newtonian hydrodynamics equations with third-order ac- 
curacy in space and second-order accuracy in time (see, e.g., 
iKifonidis etaT1l2003l and references therein). General relativis- 
tic (GR) gravity is approximated by using an "effective rela- 
tivistic potential" (for details, see Sect. 11.31 1. This approach ac- 
counts for the deeper gravitational well in GR, but still works 
in Minkowski spacetime and ignores the effects of relativistic 
kinematics. This was shown to yield excellent agreement with 
full GR calculations during the pre-explosion phas e after core 
bounce dLiebendorfer et alJ2005UMarek et~alll2006l) and will be 
demonstrated to also be a good approximation to full GR solu- 
tions of neutrino-driven winds in Sect. [2] 

The equation of state used in the simulations presented here 
is valid below densities of roughly 10 13 gcirT 3 where non- 
ideal effects due to strong interactions between nucleons can 
be safely ignored. It was used before in th e ca lculations by 
Jank a & Mullen 1 dl996h . lKifonidis et all (l2003h . and lScheck etall 
(2006). Neutrons, protons, a-particles and a representative heavy 
nucleus of the iron group (chosen to be 54 Mn) are assumed to 
be nonrelativistic Boltzmann-gases in nuclear statistical equi- 
librium. Electrons and positrons are treated as Fermi-gases of 
arbitrary degeneracy and arbitrary degree of relativity, and pho- 
ton contributions are included as well. Pressure and energy are 
corrected for Coulomb effects due to the electromagnetic inter- 
actions between nucleons and the surrounding sea of charged 
leptons. 

The transport of neutrinos and antineutrinos of all flavors is 
based on a computationally very efficient, analytic integration 
along characteristics of the frequency-integrated zeroth-order 
moment equations of the Boltzmann equation for n eutrino num- 
ber a nd energy (for details, see the Appendix of IScheck et alj 
2006). The integration yields the neutrino number and energy 
fluxes as functions of time and radius. Our approach thus ac- 
counts for the luminosity contributions due to the accretion on 
the forming neutron star. The neutrino spectra are assumed to 
have Fermi-Dirac shape with a spectral temperature that is deter- 
mined from the ratio of neutrino-energy to neutrino-number flux. 
Therefore in general the spectral temperature is different from 



the local gas temperature. The closure of the neutrino number 
and energy equations is achieved by employing the flux factor 
f(r, f) = F/(Ec), which couples the local energy (or number) 
flux with the neutrino energy (or number) density. For f(r, t) 
we use a prescribed function which was determined by fits to 
Monte Carlo transport results dJankal 19911) . This yields a reason- 
ably good approximation in the transparent and semi-transparent 
regimes but is not designed to accurately reproduce the diffusion 
limit at very high optical depths (where due to numerical reasons 
the applicability of the approach is anyway strongly constrained 
by the need of very fine grid zoning). The neutrino source terms 
in the transport equations and therefore the source terms for lep- 
ton number, energy, and momentum in the hydrodynamics equa- 
tions include the m ost relevant neutrino-matter interactions (cf. 
IScheck et al.ll2006l) . 

We note that the steepening density gradient near the neutron 
star surface requires extremely fine grid zoning for getting con- 
verged results of the neutrino-driven outflow. We typically use 
about 1000 non-equidistant radial mesh points. 

1.2. Boundary treatment 

In our simulations we replace the inner core of the neutron 
star (usually roughly 1 M e of baryonic matter) by an inner 
Lagrangian boundary of our grid, whose prescribed contraction 
is supposed to mimic the shrinking of the nascent neutron star as 
it loses energy and lepton number by neutrino emission. Using 
this inner boundary, which typically is located at a v e optical 
depth of more than 100 and a density of pib ^ 10 13 g cm 3 , does 
not only allow us to apply the simple neutrino transport approxi- 
mation described above, but also gives us the freedom to vary the 
time-evolution of the neutron star radius and of the core neutrino 
fluxes imposed at the inner grid boundary. This makes sense be- 
cause both the equation of state of hot neutron star matter and 
the neutrino transport in nascent neutron stars are not finally un- 
derstood. Changing the inner boundary conditions thus allows 
us to investigate the differences resulting from different explo- 
sion energies and timescales and from a different evolution of 
the neutrino-wind power in a given progenitor. 

Three parameters serve us to describe the motion of the in- 
ner boundary: Rf , and fo. The initial radius Ri is the radius of 
the inner core that we chose to excise from the postbounce mod- 
els we start our simulations from, Rf is the final radius of this 
core for time f — > oo, and fo is the timescale of an exponential 
contraction according to the expression 

R ib (t) =R f + 0R 1 -fl f )e-' / '° . (1) 

Our standard choice of fo = 0.1 s reproduces the contraction of 
the excised core during the first few hundred milliseconds af- 
ter bounce as foun d in full-scale supernova s imulations with the 
equation of state of lLattimer & Swestvl dl991l) . using the energy- 
depen dent neutrino trans port of the Vertex code (cf. lBuras et al.l 
2006a and also Fig. 1 in lScheck et"alll2006l) . 

In the simulations presented here we also explore the con- 
sequences of a different time-dependence of the neutrino lumi- 
nosities imposed at the inner grid boundary (see Sects. [XT] and 
|3.5t . The explosion energy of a model is mostly determined by 
the choice of the initial values of these luminosities (in particu- 
lar those of v e and v e ). These initial values are constrained by the 
prescribed total loss of neutrino energy from the core during the 
proto-neutron star cooling, AE 1 "^, and by t he total loss of lep- 
ton number, AYe iCore (see Eqs. D.13-D.16 in lScheck et al.| [2006, 
where Ly Ot ' fL in the last equation should be replaced by AE™^). 
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The relative contribution of v e to the total core luminosity is set 
to 20%, (i.e., K v , = 0-2 in terms of the parameters introduced 
in IScheck etal]l20 06). the contribution of v e is determined from 
requesting AF e coie = 0.3, and the muon a nd tau neutrino con - 
tributions then follow from Eq. (D.12) in IScheck et all (|2006). 
The mean energies of the neutrinos entering the computational 
grid at the inner boundary are chosen to be (e Vc ) lb = 12MeV, 
(e p ) ib = 16MeV, and {e Y j h = 20MeV when v x denotes muon 
and tau neutrinos and antineutrinos. These energies are kept con- 
stant with time. 

Because of the contraction and postbounce accretion of the 
proto-neutron star, the density and optical depth in the layers 
near the inner grid boundary can increase to such large values 
that the application of our transport approximation becomes in- 
efficient by the required very fine zoning, and the equation of 
state fails to describe the dense stellar matter. Whenever the v e 
optical depth begins to exceed a certain value (usually chosen to 
be 300), we shift the inner boundary to a larger radius /?ib(? C ut) 
and thus to a larger mass shell where the neutrino optical depth 
is significantly lower (usually 200 for v e ). The additional excised 
baryonic mass is added to the previous core mass and the gravi- 
tational mass of the new, increased core is set equal to the grav- 
itational mass computed at radius /?ib(? C ut) where the new inner 
grid boundary is placed (see Sect. 11.3b . The subsequent motion 
of the new boundary for t > f cut is assumed to follow the function 

=Rf+ (tfibftat) - «f) exp [v(t - t cut )/(R ib (t cm ) - R f )] , (2) 

where v < is the recession velocity of the mass shell of the 
new boundary at time f cut . The new boundary contracts in a very 
similar way as the previous one because the removed shell is 
very narrow. The neutrino luminosities and mean energies of 
the streaming neutrinos imposed at the new boundary at t — f cut 
are chosen to be the values computed with the transport scheme 
at this radius and to have the same time behavior as the initial 
boundary luminosities and mean energies. 



1.3. Gravity 

Relativistic effects are taken into account in our Newtonian hy- 
drodynami cs code by using an "ef fective relativistic gravitational 
potential" (Ra mpp & Jankall2002h . The simulations pre sented in 
this paper (different from those of IScheck et al.ll2006l) employ 
the im proved version of this potential described by iMarek et al.l 
(l2006h . who found excellent agreement with fully relativistic 
calculations during core collapse and the first several hundred 
milliseconds after core bounce (tests for the later neutrino-wind 
phase can be found inSect.[2]l. 

Ac cording to iRampp & Jankal (|2002) and IMarek et all 
(2006), the relativistic equation of motion can be rearranged in a 
form similar to the Newtonian one by replacing the Newtonian 
gravitational potential with a modified TOV potential: 



with the metric function 



•DtovM = -AnG 



Cdr 
J ~' 



dr' I m TO v r' 3 P\ 1 / pc 2 + e + P 



\ An 



r 2 



pc^ 



(3) 



where p is the rest-mass density, e — pe the internal energy den- 
sity with e being the specific internal energy, and P the gas pres- 
sure. The usually rather small corrections of the gravitational 
poten tial due to neutrin o pressure, energy density, and flux terms 
(see lMarek et al.ll2006l) are neglected in Eq. 0). The "modified 
TOV mass" m-rov is given by 



m TO v(/) = An J dr'r' 2 (p + -^jr 



1 + - 



2Gmjoy 



(5) 



The extra factor F in Eq. (0]l compared to the relativistic defi- 
nition of the TOV mass enters the mass integral for reasons of 
consistency with the Newtonian hydrodynamics equations and 
accounts for the fact that in the Newtonian code there is no dis- 
tinction between local proper volu me and coordinate volume, 
i.e., d'V = dV (for more details, see lMarek et alj|2006l) . 

There is, however, an important diff erence of ou r calcu - 
lations compared to those performed by Mar ek etaTl ([2006 ). 
While the latter included the whole neutron star down to the cen- 
ter, the use of the inner grid boundary at a radius Rn, > in the 
present work prevents the evaluation of the integral in Eq. (0| 
within the neutron star core. We solve this problem by starting 
our calculations with a given value of the modified TOV mass 
of the core at t = 0, m T ov(^ib> 0), which was provided to us as 
part of the data set for the initial conditions of our simulations. 
For t > we then approximately evolve the modified TOV mass 
according to the expression 



m TO v(^ib, = m TO v(Rib, 0) 



f 4V)d/' 
Jo 



-X 

Jo 



AnR 2 b {t')P ih (t')^dt' , (6) 



where the second term on the rhs yields the energy loss from 
the neutron star core by the total neutrino luminosity at the inner 
boundary, L^if), and the last term represents the compression 
(PdV) work done on the core at the contracting inner boundary. 
The total modified TOV-mass at radius r, which we consider as 
"gravitational mass", is thus given by 



rnxovCn = mjov(R ib , t)+An \ dr'r' 2 [p + — ) F . 

Jft,b v c ' 



(7) 



(4) 



2. Comparison with fully relativistic wind solutions 

In this section we will discuss our simulation approach in view 
of other published work on relativistic steady-state solutions 
for neutrino-driven winds. General relativistic (GR) effects have 
been recognized to cause important changes of the neutrino wind 
properties, e.g., to lead to a decrease of the ex pansion timescale 
and to an increase of th e wind e ntropy, see Qian & Wo oslevl 
(fl996h . ICardall & Fuller! (fl997h . ISumivoshi et alJ (120001) . A 
comprehensive discussion of these effects in comparison with 
the Newtonian treatment w as provided by lOtsuki et al.l ([2000) 
and lThompson et al.l d2001l) . 

In our simulations we account for relativistic gravity only by 
using a modified effective potential (Sect. fOl ), but otherwise we 
solve the Newtonian equations of hydrodynamics. We also ig- 
nore relativistic redshift and ray bending effects in o ur descrip- 
tion of the neutrino transport (cf. lScheck et alj|2006l) . 

The use of the generalized potential in a Newtonian hydro- 
dynamics code was shown previously to yield results in very 
good agreement with relativistic core -collapse simulations up 
to several 100 ms after core bounce ( Liebendo rfer et all [2005; 
Mare k et ai1l2006l) . For the later neutri no-wind phase we tried 
to compare with solutions plotted by Thomps on et al.l (2001) 
for cases when our neutron star masses, neutron star radii, and 
neutrino-heating rates were similar to the ones considered in that 
paper. Unfortunately, we were unable to find moments in our 
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Fig. 1. Derivatives of the velocity (top), density (middle), and 
temperature (bottom) as functions of radius from our hydrody- 
namic model M15-ll-rl at 1.5 s after bounce (black lines) com- 
pared to these derivatives as computed from the relativis t ic sta- 
tionary wind equations of Eqs. (5)-(7) in lThompson et al.ld200lh 
(blue curves). These equations were evaluated by using the val- 
ues of all gas quantities as provided by our hydrodynamic model. 
Consistency between our hydrodynamics results (with approxi- 
mative treatment of relativity) and the fully relativistic wind so- 
lution would require the corresponding lines to lie on top of each 
other. The agreement is very good and the sonic point is located 
at about 80 km in both cases. This location is a critical point of 
the wind equations, which explains the pathological behavior of 
the curves there. The evaluation of the expression for the tem- 
perature gradient is numerically inaccurate in a region where the 
two terms in Eq. (7) of Thompson et al.'s paper are very large 
and have opposite signs, in which case dT/dr becomes slightly 
positive while the hydrodynamical result is still negative. 



simulations where al l relevant parameters m atch up exactly the 
cases considered by Thomp son et al.l (1200 ll) . As far as a com- 
parison was possible, we observed satisfactory agreement in the 
main properties characterizing the wind. 

A more quantitative comparison is hampered by the fact that 
relativistic neutrino- wind simulations are not available to us. We 
therefore decided to make use of Eqs. (5)-(7) for the velocity 
derivative, dv/dr, the density derivative, dp/dr, and the temper- 
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Radial profiles of the net total neutrino heating rate q tot 
\/ e and v e luminosities L Y (second panel, in bethe per sec- 



Fig. 2. 

(top), 

ond or 10 51 ergs : ), mean neutrino energies (e v ) (third panel; 
Eq. HJ, rms energies (e v ) rms (fourth panel; Eq. |9|, and flux fac- 
tor for model M15-ll-rl at 1.5 s after bounce. For comparison 
with Thompson at al. (2001), we also show (e v ) and (e v ) lms as 
used in that paper, rescaled to our values of L Vc and Ly c , the flux 
factor with and without relativistic corrections from that work, 
and the Newtonian and GR charged-current heating plus cool- 
ing rates using Thompson et al.'s formulas, evaluated with our 
neutrino luminosities and rms energies (red and blue curves) or 
with the rescaled rms energies of Thompson et al. (green curve). 
The vertical dashed line marks the position of the v e -sphere, and 
the vertical dotted line the neutron star "surface" at a density of 
10" gem" 3 . 
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ature derivative, dT /dr, in the paper of Thomp son et al.l (|2001). 
Figure Q] shows these derivatives as functions of radius at a cer- 
tain time for one of our models, compared to the results from 
Thompson et al.'s fully relativistic expressions. Evaluating the 
latter, we took all quantities on the rhs. of the formulas (veloc- 
ity v, adiabatic sound speed c s , density p, gravitational mass M, 
neutrino heating rate q, etc.) from our model. Ideally, the pairs of 
corresponding curves in Fig. \T\ should fall on top of each other, 
which would demonstrate consistency of both calculations. The 
overall agreement of the two cases is very good, with a small 
difference being visible only around the maximum of the ac- 
celeration, which, however, is located at the same radius. Also 
the sonic point is nearly at the same position of about 80 km 
(we are not disturbed by the pathological behavior of the curves 
in this region, where the expressions for the derivatives have a 
critical point). We therefore conclude that our approach repro- 
duces the most important features of the relativistic solution, and 
that relativistic kinematics (which we ignore) is of minor impor- 
tance compared to the effects of the stronger GR potential, which 
makes the proto-neutron star more compact and the density and 
temperature gradients in the neutrinospheric region steeper than 
in Newtonian gravity. 

We also compar ed our neutri n o hea ting a nd cooling rates 
with t hose used by lOtsuki et"aLl d2000l) and IThompson et al.l 
(2001). Figure [2] shows the radius-dependent net (i.e., heating 
minus cooling) specific rate of neutrino energy deposition by the 
^-proc esses according to Eqs. (20) and (21) of IThompson et all 
(2001) with and without corrections for relativistic redshift and 
ray bending, evaluated at all radii with the stellar parameters 
and the neutrinospheric values of the v e and v e luminosities and 
mean energies from one of our simulations. The data were taken 
from the same model and time used in Fig. Q] The behavior 
of bot h curves agrees qualitatively with Fig. 5a of Otsuki ^tail 
(2000). Close to the neutrinosphere ray bending effects enhance 
the net heating (since GR causes a reduction of the flux factor 
as visible for Thompson et al.'s prescription of this quantity in 
Fig. |2), whereas gravitational redshifting of the neutrino lumi- 
nosities and energies grows monotonically with distance from 
the neutrinosphere an d finally wins, red ucing the GR rate below 
the Newtonian value. lOtsuki et all d2000l) performed test calcu- 
lations to disentangle the influence of GR corrections in the neu- 
trino treatment from that of the relativistic terms in the wind 
structure eq uations. In spite of the sizable change of the local 
heating rate. lOtsuki et al.l d2000h found that neutrino redshift and 
ray bending have only little im pact on the wind entrop y. Similar 
conclusions were arrived at by Thomp son et al.l d200ll) . 

In Fig. |2] also the total specific rate of neutrino heating and 
cooling from our hydrodynamical model is displayed. This rate 
includes all contributing processes, i.e. besides the /J-reactions 
of v e and v e absorption and production also energy transfer by 
the scattering off electrons, positrons, and free nucleons, and 
neutrino-antineutrino pair annihilation, to which neutrinos of all 
flavors contribute (cf. the appendix of Sche ck et al.l l2006). This 
total rate is similar to th e neutrino captur e and emission rates 
of Eqs. (20) and (21) of Thomp son et al.l d200ll) . because for 
the considered situation the neutrino luminosities are high (see 
Fig. |2) and therefore the wind mass loss rate is large and the wind 
entropy fairly low (s w 55 50 k% per nucleon). At such conditions 
of high wind density and modest abundance of e + e~-pairs, the 
other reactions do not contribute significantly to the total rate of 
energy deposition. 

We point out here that our approximative treatment of neu- 
trino transport evolves the transport solution self-consistently 
with the temperature and density structure of the stellar medium. 



This is different from the light-bulb approach of previous steady- 
state or hydrodynamical wind studies (e.g.. ISumivoshi etalj 
2000: lOtsuki et al.ll2000t IThompson et alj200ll) . Inside the neu- 
trinosphere neutrinos and matter are in equilibrium, around the 
neutrinosphere neutrinos begin to decouple thermodynamically 
from the medium, and at some larger distance they start stream- 
ing freely. A changing radial structure of the contracting neutron 
star leads to changes of the neutrino luminosities and mean en- 
ergies, and the gradual loss of neutrinos drives the cooling and 
deleptonization of the surface-near layers of the neutron star. 
In previous wind studies (except full supernova models), such 
a coupling and interdependence was ignored. Close to the neu- 
tron star surface the flux factor (or flux dilution factor) used in 
our transport , which is based on a Monte Carlo calibration by 
Janka (1991), is lower tha n the vacuum app roximation chosen 
bv IThompson et al.l (1200 ll) and lOtsuki et al.1 (120001) , see Fig. |2] 
Thompso n et al.l (12001 1) have tested the improved description by 
iJankal (Il99ll) and found that its effects are negligible for the 
range of model conditions considered by them. This, however, 
is true only during phases where the density gradient near the 
neutrinosphere is very steep and in regions where the neutrino 
luminosities have already reached their asymptotic values. 

The most important difference of our simulations compared 
to other relativistic wind studies is the different treatment of 
the spectra in the neutrino transport. In our "grey" but non- 
equilibrium description of neutrino number and energy trans- 
port, we determine a neutrino spectral temperature that is inde- 
pendent of the matter temperatur e and can be diffe rent from it 
(for details, see the appendix of IScheck et al.l 120061) . This leads 
to higher mean energies of v e and v e radiated from the neutri- 
nosphere than considered in the other works. Figure |2] shows 
these mean energies as functions of radius, defined once as the 
ratio of neutrino energy flux to neutrino number flux, 



(e v ) = 



and another time as rms energy, 



( f v)rms — ksTy 



T5(T] V ) 



(8) 



(9) 



which is the energy which enters the calculatio n of the neu- 
trino absorption rates on nucleons (cf., for example. lScheck et al.l 
2006). In Eq. (|9), T Y and j] v are the spectral temperature and de- 
generacy, assuming that the neutrino spectra have Fermi-Dirac 
shape, in which case 



r„(y) 



-r 

Jo 



dx ■ 



x" 



1 + exp (x - y) 



(10) 



For comparison, Fig. |2] also presents the corresponding mea n 
energies and rms energies as used bv iThompson et all d200ll) . 
appropriately scaled by l)J A to account for the larger neutrino 
luminosities considered here, and taking tj v — for neutrinos 
and antineutrinos. The net heating rate computed with these rms 
energies is significantly lower than the heating rate from our hy- 
drodynamic model (Fig. [2). 

At first glance, our mean energies for v e and i> e might appear 
on the large side. One must, however, take into account that the 
mean energies in our simulations are significantly lower in the 
first ~ 0.5 seconds when the neutron star is still rather extended, 
and only increase as it heats up during contraction. They reach 
a maximum between one and two seconds after bounce to de- 
crease afterwards as the proto-neutron star cools (see Sect. [3). 
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Table 1. Model parameters of our spherically symmetric simula- 
tions. The different models are characterized by the chosen con- 
traction of the inner grid boundary, which is expressed in terms 
of the final radius Rf and the exponential contraction timescale 
fo (cf. Eq. [TJ. Different choices of these values are indiced by 
the extensions "rl", "r2", etc. of the model names. In addition, 
different initial luminosities of v e plus v e (measured in bethe [B] 
= 10 51 erg per second) are imposed at the inner grid boundary in 
case of our standard luminosity behavior (constant until 1 s and 
then a f~ 3 ^ 2 decay). These variations are reflected by the exten- 
sions "11", "12", etc. in the model names. Moreover, the time- 
dependence of the boundary luminosity has been modified to a 
luminosity decay that is more rapid than in the standard descrip- 
tion (Eqs.QT|and[T2l models with "It" in their names). 



Model 


Contraction 
(Rf, k) 


[B/s] 


Progenitor Mass 
[M ] 


M15-ll-rl 


9km; 0.1s 


52.5 


15 


M15-ll-r2 


9 km; 0.2 s 


52.5 


15 


M15-ll-r5 


11km; 0.1s 


52.5 


15 


M15-ll-r6 


14 km; 0.1s 


52.5 


15 


M15-12-rl 


9km; 0.1 s 


38.6 


15 


M15-13-r3 


10km; 0.1s 


35.8 


15 


M15-lt2-r3 


10km; 0.1s 


55.2 


15 


M15-ltl-r4 


10.5 km; 0.1s 


55.8 


15 


M10-ll-rl 


9km; 0.1 s 


52.5 


10 


M10-15-r3 


10km; 0.1s 


30.3 


10 


M20-ll-rl 


9km; 0.1 s 


52.5 


20 


M20-13-r3 


10km; 0.1s 


35.8 


20 


M20-14-r3 


10km; 0.1s 


33.1 


20 


M25-15-r4 


10.5 km; 0.1s 


30.3 


25 



So the conditions plotted in Fig. [2] correspond to a time when 
the neutrino luminosities are still rather high and the mean en- 
ergies in this phase at their maximum. Moreover, one should re- 
member that we ignore gravitational redshifting in our transport. 
The redshift from the neutrinosphere at radius R v to infinity after 
the contraction of the neutron star can become quite significant, 

Vl - 2GM/(R v c 2 ) = yJl-R s /R v * 0.7. ..0.8 for \ > RJR V > ±, 
which reduces the mean energies for a distant observer by typ- 
ically 20-30%. In our Newtonian transport treatment we prefer 
to use the higher neutrinospheric energies for evaluating the neu- 
trino heating, because the neutrino-wind properties are mostly 
determined by the heating just outside of the neutrinosphere, 
where it is also strongest. 



3. Results 

In this section we will present the results of our time-dependent 
hydrodynamic simulations. An overview of the computed mod- 
els will be given in Sect. 13.11 We will begin with the descrip- 
tion of a reference case in Sect. 13.21 This model was computed 
for a certain choice of the time-dependent contraction of the in- 
ner grid boundary and of the neutrino luminosities imposed at 
this boundary. We will then demonstrate that basic features of 
the wind termination shock can be understood by simple ana- 
lytic considerations (Sect. |3.3l l. Next, we will discuss the neu- 
trino wind evolution in different progenitor stars (Sect. 13.41 ), and 
finally will investigate the influence of variations of the condi- 
tions at the inner boundary (Sect. [331 ). 



3.1. The computed models 

A list of computed ID models with their characterizing parame- 
ters is given in Table [TJ We have performed simulations for pro- 
genitor stars with 10.2 M (data provided b y A. Heger, personal 
comm unication), 15 M Q (m odel s!5s7b2; Wooslev & Weaver] 
119951) . 20 M (model s20.0, 



Wooslev et al. 



2002]), and 25 M 



(model s25a28. iHeger et aT]|2001l) . Pre-collapse pr ofiles of the 
last thr ee models are plotted in the appendix of iBuras et ail 
(l2006al) . The progenitor cores were computed through core col- 
lapse and bounce with the neutrino-hydrodynamics code Vertex 
and provided to us as initial conditions for the present studies 
a few milliseconds after shock formation (A. Marek, personal 
communication). Extensions of the model names ("rl", "r2",....) 
indicate different prescriptions for the contraction of the inner 
grid boundary, whose motion was varied by choosing different 
values of the final radius Rf and of the exponential contraction 
timescale fo (Eq. [TJ. Larger numbers in this sequence corre- 
spond to less quickly contracting or less compact neutron stars. 
Moreover, we varied the sum of the v e and v e luminosities im- 
posed at the grid boundary with respect to the initial value as 
well as time-dependence. In most of the calculations the lumi- 
nosities were chosen to be constant during the first second of 
postbounce evolution and to de cay proportional to r 3 ^ 2 after- 
wards as in Schec k~et al.l d2006l) . Such models are labelled by 
the extensions "11", "12", etc., with a higher number meaning a 
lower initial value of the boundary luminosity. In another set of 
calculations the boundary luminosities were assumed to have a 
smoother time-dependence (with no jumps in the time deriva- 
tive) and in particular with less neutrino energy radiated at late 
postbounce times. The luminosities were prescribed as 



L(t) 



with 



/(?) 



(La 

Uo. 



if t < 0.5 s ; 
,f{t) , if t > 0.5 s , 



exp [-it - 0.5) 2 ] + b [ 1 + (t - 0.5)" ]" 
d+b) 



(ID 



(12) 



where the time t is measured in seconds. The corresponding 
models can be recognized by the letters "It" in their names. The 
parameter n was set to 1.5 in both cases, while b - 0.2 was used 
for model M15-ltl-r4 and b = 0.3 for M15-lt2-r3. 

A comparison of these models allows us to study the influ- 
ence of different contraction behavior of the nascent neutron star. 
The contraction determines the release of gravitational energy 
from the mantle layers of the compact remnant. The accretion 
luminosity generated in the mantle adds to the core flux (given 
by the imposed boundary condition) and has an influence on the 
explosion timescale and explosion energy of a model and thus on 
the location of the mass cut and the baryonic mass of the neu- 
tron star. The corresponding gravitational mass, which decreases 
when energy is lost in neutrinos (Eq. [6j, the radius of the neu- 
tron star, and the luminosities and mean energies of the radiated 
neutrinos are crucial parameters that dir ectly affect the neutrino - 
wind properties as functions of time (see Oi an & Wooslevll 996). 

We note that the supernova models we study here do not per- 
mit us to change individually and independently all parameters 
and conditions that affect the neutrino wind properties and that 
determine the behavior of the wind termination shock. The wind 
depends, e.g., on the neutron star gravitational potential and thus 
on the neutron star mass. The latter becomes larger when the 
postbounce accretion phase lasts longer and the explosion hap- 
pens later, or when the progenitor is more massive and therefore 



8 



Arcones et al.: Neutrino-driven supernova outflows 



the iron core and postbounce accretion rate are larger. More mas- 
sive progenitors thus tend to produce neutron stars with bigger 
masses. For this reason one cannot disentangle the influence of 
the progenitor structure on the wind termination shock from the 
effects of the neutron star mass on the neutrino wind. 

In order to structure the discussion, we decided to first de- 
scribe basic features in case of a 15M Q reference model, then 
to vary the boundary conditions for this model, and finally to 
present the results for different progenitors. 

3.2. Evolution of a reference case 

In our reference model, M15-ll-rl, the conditions at the inner 
boundary were chosen such that the 15M star explodes with 
an energy of ~1.5xl0 51 erg and the neutron star attains a final 
gravitational mass of ~ 1.2 M (1.4 M Q baryonic) and a radius of 
10 km (Tabled). 

A mass-shell plot for the space-time evolution of this model 
is given in Fig. [3] The explosion sets in about 200 ms after 
bounce (at the time f exp given in Table [2] which is defined as 
the moment when the total energy of expanding matter starts to 
exceed 10 49 erg). At this time the stalled shock is revived by neu- 
trino heating and starts continuous expansion with an average 
velocity of roughly lO.OOOkms -1 . On its way out the shock re- 
verses the infall of the swept-up matter. After the onset of the ex- 
plosion, ongoing neutrino energy transfer drives an outward ac- 
celeration of heated material in the gain layer around the neutron 
star. At the interface between this dilute neutrino-driven wind 
and the denser outer ejecta a contact discontinuity is formed. 
Even farther behind the forward shock, the neutrino-driven wind, 
whose velocity increases rapidly with distance from the neutron 
star, collides with more slowly moving material and is deceler- 
ated again. The strongly negative velocity gradient at this loca- 
tion steepens into a reverse shoc k when the wind veloc ity begins 
to exceed the local sound speed dJanka & Miillerl 1995b . First in- 
dications of a forming wind termination shock can be seen in 
Fig.[3]at t 350 ms post bounce at a radius r ~ 300 km. 

Figure |4] displays the v e and v e luminosities and the mean 
energies emitted by the nascent neutron star (gravitational red- 
shift effects are ignored). One can see the accretion phase with 
its production of accretion luminosity ending at the time the ex- 
plosion sets in. The following plateau phase until t » 1 s and 
subsequent decay of the luminosities show the influence of the 
time-dependence of the imposed boundary fluxes. This is also 
the case for the mean neutrino energies. Their values increase 
during the first second of postbounce evolution because the in- 
ner grid boundary and the neutron star radius contract (Fig. [3]). 
Consequently, the outer layers of the neutron star heat up due to 
the conversion of gravitational energy to internal energy by com- 
pression. After one second the rapid contraction is over and the 
decay of the boundary luminosities leads to less energy transport 
into these layers, which therefore begin to cool down, causing 
the mean energies of the radiated neutrinos to decline. 

Figure [4] also provides information about the total energy 
carried away by neutrinos and antineutrinos of all three lepton 
flavors, AEtot, as a function of time, and the corresponding re- 
duction of the gravitational mass of the nascent neutron star. The 
latter is taken to be the modified TOV mass (Eq.|7]) at the neutron 
star radius R ns (which is defined as the radius where the density 
is 10 n gcnT 3 ). In contrast, the baryonic mass of the neutron 
star, given by the rest mass enclosed by the radius R ns , initially 
increases in the course of accretion. After the explosion has been 
lauched, it decreases again only slightly due to the mass loss in 
the neutrino-driven wind. 



The middle panel on the right side of Fig.|4]reveals that only 
about 50% of the explosion energy (defined as the integral of 
the total internal, kinetic, and gravitational energies in all zones 
where the sum of these energies is positive) are carried by the 
neutrino-heated shell of matter expanding right behind the shock 
after the onset of the explosion. The rest is contributed by the 
early neutrino wind, and after 2 s the energy has reached 95% of 
its final value. The panel below gives the expansion timescales 
of the ejected mass shells in the neut rino wind. The first defi- 
nition follows Oian & Woosley (1996, Eq. 60), who introduced 
the dynamical timescale as 



^"dyn — 



fc,r=0.5MeV 



(13) 



We co mpare this with the cooling timescale used by Otsuk{eTaTj 
(120001 Eq. 23), 



f 

J k: 



■jt B r=0.5MeV/e 



dr 



(14) 



B r=0.5MeV 



variations of which were considered bv lWitti et al.l d 19941 cool- 
ing time between T = 7 x 10 9 K and T = 3 X 10 9 K) and 
WanajoetaL|{2001, cooling time between k^T — 0.5 MeV and 
k^ T = 0.2 Me V). The th ird definition we consider is the one 
of Thompson et al. (2001, Eq. 32), who employed the e-folding 
time of the density instead of that of the temperature, 



Idp 
p dr 



(15) 



i B r=0.5MeV 



where we set for our Newtonian simulations y — 1 in Thompson 
et al.'s Eq. (32). As can be expected from the fact that the 
wind is radiation-dominated and therefore s oc T 3 /p ~ const, 
the timescale r p is always significantly shorter than the cooling 
timescale tj (Fig.|Ul. Ideally, in such a situation one would ex- 
pect Tt/t p = 3, which is more closely reached at later stages 
when the wind entropy is higher (Fig. |5). Due to the different 
mathematical expressions in Eqs. (fT4l and (fT~5b . the factor 3 is 
never exactly realized. The third timescale, Eq. ( fT3l , yields a re- 
sult that is between the other two values during most of the com- 
puted postbounce evolution and comes closer to the timescale of 
Eq. ( fT4b in the late stages of the simulations. 

In Fig. [5] the radial profiles of different wind quantities are 
given for our reference model M15-ll-rl at a number of post- 
bounce times. The neutrino heating accelerates the wind to a 
peak velocity of about 25% of the speed of light for neutrino 
luminosities L Vc ~ Ly c « 3 x 10 52 ergs~'. The maximum ve- 
locity decreases as does the heating rate when the luminosities 
decline with time. The density and temperature in the wind re- 
gion follow roughly the usual r~ 3 and r _I behavior, respectively, 
in the region where the entropy is a constant. The profiles are 
slightly steeper and thus closer to these power laws at later times 
when the wind entropy is higher and the wind therefore more 
dominated by radiation pressure. The radial profiles of p and T 
also steepen at larger distance from the neutron star, leading to 
a visible increase of the wind acceleration at the point where 
free nucleons recombine to a-particles and the neutrino heating 
ceases (see the corresponding panels in Fig. [5]). At this radius 
the entropy of the outflow reaches its final value. This asymp- 
totic wind entropy increases from about 60 &b per nucleon at 1 s 
to around 90 k% at 10 s. 

For (approximately) the same values of the v e and v e lumi- 
nosities, model M15-ll-rl tends to yield somewhat lower ex- 
pansion timescales, slightly lower entropies, and a bit higher 
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Fig. 3. Mass shell plot for the evolution of 
model M15-ll-rl. The explosion occurs about 
0.2 seconds after bounce. The red line marks 
the supernova shock, the blue line the wind 
termination shock, the orange lines the loca- 
tions where the expansion velocity of the gas 
equals the local sound speed (sonic points), 
and the green line the neutron star radius de- 
fined as the location where the density drops 
below 10" gcirr 3 . A contact discontinuity sep- 
arates the dense shell of ejecta that were accel- 
erated by the outgoing shock from the very di- 
lute neutrino-driven wind. Mass shells in the 
wind are labelled by the corresponding en- 
closed baryonic masses. 




Fig. 4. Luminosities, mean energies according to Eq. ([8j, and rms energies (Eq.[9] (e) rras > (e)) of v e and v e , and total energy radiated 
in neutrinos of all flavors for model M15-ll-rl as functions of time (left), measured outside of the nascent neutron star (at a radius 
of 500 km). Note that we do not include gravitational redshifting in our neutrino treatment. The rapid decline of the luminosities 
after about 0.2 marks the end of the accretion phase of the forming neutron star at the onset of the explosion. The panels on the rhs 
side give the baryonic mass and the gravitational mass (Eq. |7J of the neutron star in model M15-ll-rl, the explosion energy, and 
the expansion timescales of the neut rino-driven wind as func t ions of time. For the latte r, the results from three different definitions 
are displayed, namely those used bv lOian & Wooslevl (119961) . lOtsuki et at] (12000). and Tho mpson et al.l yOOl), given in Eqs. ( TOl , 
( TBI , and ( 115) , respectively. 

mass loss rates than those found by Thomps on et al.l d2001l) . see the one hand by the smaller gravitational mass of the neutron star 
Figs. 5, 8, and Tables 1 and 2 there. This can be understood on in our model compared to the canonical 1.4 Mq star considered 
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Fig. 5. Radial profiles of the net neutrino-heating rate q, density, temperature, entropy (left, from top to bottom), velocity, mass 
loss rate, electron fraction Y e , and mass fractions of free protons and a particles for the neutrino wind in model M15-ll-rl at 
different postbounce times. For the mass fractions only the information for the first and last moments of time is provided. The wind 
termination shock is clearly visible in its effects on the velocity, density, temperature, and entropy of the outflow. 



bvlThompson et al. (2001), and on the other hand it is caused by 
our larger heating rates due to the higher mean neutrino energies 
(cf. our discussion in Sect. [2]). These differences affect the char- 
ac teristic wind parameter s with different sensitivity. According 
to lOian & Wooslevl dl996l) . the entropy scales with the neutrino 
luminosity L, the mean neutrino energy e, the neutron star radius 
R, and the neutron star mass M like 

s oc L- 1,6 <e- 1/3 R- 2/3 M , (16) 
the expansion timescale like 

rccL^e^RM, (17) 
and the wind mass loss rate like 

M oc L 5/ V 0/3 /? 5/3 AT 2 (18) 

(modificatio ns of these relations du e to relativistic effects were 
addressed bv lThompson et al1l2001l) . 

During the first ~ 2 seconds after the onset of the explo- 
sion, the neutrino wind is p-rich, i.e. 5^ > 0.5. This is in 



agreement with explosion models that e mploy a Boltzmann 
solver for the sp ectral neutrino transport (see Bu ras et al.l l2006b; 
iPruet et ail 120051) . Afterwards the electron fraction drops be- 
low 0.5, and gradually the wind develops increasing neutron 
excess. Qualitatively, this trend to lower Y e at later times is 
reproduced when the neutrino luminosities and mean energies 
from the simulation are inserted into the simple analytic rela- 
tion Y e ~ [ 1 + (L Vt e v JL Vc e Ve )] ~ l , although the values do not 
agree quantitatively. We emphasize here that the gray and ap- 
proximative treatment of the neutrin o transport emplo yed in this 
work (for a critical assessment, see IScheck et al.ll2006l) is also 
not able to yield reliable results for the electron fraction in terms 
of absolute numbers. The competition of v e and v e absorption on 
free neutrons and protons sensitively determines the asymptotic 
value of Y e , an accurate calculation of which requires detailed in- 
formation of the neutrino and antineutrino spectra in the comov- 
ing frame of the expanding wind matter. The wind at late times 
might therefore become significantly more neutron rich than pre- 
dicted in our models. Figure [5] also reveals that the mass loss 
rate reaches its asymptotic value closest to the neutrinosphere, 
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Fig. 6. Radius, density, entropy (left, from top to bottom), velocity, temperature, and electron fraction Y e as functions of time along 
the trajectories of different mass shells that are ejected in the neutrino-driven wind of model M15-ll-rl. The times correspond to 
the moments when the mass shells cross a radius of 100 km. After a very rapid expansion, the wind is abruptly decelerated by the 
termination shock. This leads to an increase of the entropy by more than a factor of two, and to a subsequently much slower decline 
of the temperature and density. 



Table 2. Results of the ID models at the end of the simulations at t = 10 s after bounce. Mb ar is the baryonic mass of the neutron 
star, M erv its gravitational mass (Eq. [7}. Both masses are computed for the matter inside the neutron star radius R ns . This radius is 
defined as the location where the density is 10 1 1 g cirT 3 . M e tf denotes an "effective mass" of the neutron star, for which a Newtonian 
force equals to the gravitational force associated with the modified TOV potential of Eq. (0 at radius R ns . AE tot is the total energy 
radiated in neutrinos of all flavors (measured in bethe [B] = 10 51 erg), L Vc and Ly c are the luminosities of electron neutrinos and 
antineutrinos measured at 500 km (without gravitational redshifting), (e v .) and (ej> o ) are the corresponding mean energies, E exp is the 
explosion energy (note that this energy can still decrease somewhat after 10 s because of the negative binding energy of the outer 
stellar layers, which is not included in the given numbers), t &xp is the postbounce time when the explosion sets in (defined as the 
moment when the energy of expanding postshock matter exceeds 10 49 erg), s w is the asymptotic wind entropy per nucleon, and s ls 
the entropy of the outflow after its deceleration in the wind termination shock. 



Model 


time 


M bar 


M s „ 


M cff 


A£ tot 












F 

^exp 


►exp 


'^wind 


Sis 




[s] 


[M B ] 


[M e ] 


[M s ] 


[100B] 


[km] 


[B/s] 


[B/s] 


[MeV] 


[MeV] 


[B] 


[s] 


[£ B /nuc] 


[£ B /nuc] 


M15-ll-rl 


10.0 


1.399 


1.207 


1.910 


3.422 


10.09 


1.73 


2.06 


16.19 


19.52 


1.486 


0.201 


91.10 


190.55 


M15-ll-r2 


10.0 


1.446 


1.263 


2.031 


3.335 


10.27 


1.96 


2.28 


16.37 


19.50 


1.174 


0.341 


94.19 


191.29 


M15-ll-r5 


10.0 


1.394 


1.208 


1.703 


3.307 


12.94 


1.78 


2.11 


14.57 


17.62 


1.371 


0.221 


71.97 


131.94 


M15-ll-r6 


10.0 


1.440 


1.258 


1.648 


3.274 


16.71 


1.61 


1.93 


13.04 


15.98 


1.043 


0.241 


61.04 


83.63 


M15-12-rl 


10.0 


1.473 


1.280 


2.116 


3.451 


9.94 


1.74 


2.12 


16.46 


20.30 


1.019 


0.381 


100.14 


193.14 


M15-13-r3 


10.0 


1.545 


1.341 


2.121 


3.709 


11.22 


1.99 


2.34 


15.88 


18.97 


0.709 


0.701 


92.75 


155.84 


M15-lt2-r3 


10.0 


1.397 


1.253 


1.906 


2.602 


11.28 


1.33 


1.54 


15.04 


18.12 


1.239 


0.221 


90.19 


132.57 


M15-ltl-r4 


10.0 


1.395 


1.260 


1.878 


2.421 


11.81 


0.96 


1.12 


14.33 


17.40 


1.231 


0.221 


91.02 


96.92 


M10-ll-rl 


10.0 


1.314 


1.132 


1.721 


3.251 


10.22 


1.63 


1.96 


15.85 


19.06 


1.247 


0.321 


83.48 


476.93 


M10-15-r3 


10.0 


1.344 


1.187 


1.745 


2.817 


11.49 


1.47 


1.72 


14.98 


18.04 


0.716 


0.421 


80.94 


353.32 


M20-ll-rl 


10.0 


1.422 


1.233 


1.973 


3.388 


10.10 


1.74 


2.05 


16.23 


19.52 


1.486 


0.181 


94.18 


127.89 


M20-13-r3 


10.0 


1.595 


1.411 


2.310 


3.383 


11.10 


1.60 


1.90 


15.66 


18.79 


0.375 


0.761 


105.10 




M20-14-r3 


10.0 


1.523 


1.332 


2.118 


3.437 


11.03 


1.61 


2.00 


15.69 


20.00 


0.847 


0.421 


95.32 


106.84 


M25-15-r4 


10.0 


1.971 


1.657 


2.944 


5.924 


11.56 


2.95 


3.58 


16.88 


20.01 


1.700 


0.401 


113.75 


117.89 
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and only then Y e and finally the entropy reach their asymptotic 
values. 

At a radius of a few 1000 km, the supersonic wind is abruptly 
decelerated down in the termination shock. The compression 
leads to a density and temperature increase. The conversion of 
kinetic to internal energy in the shock boosts the entropy to 
more than twice the wind entropy in model M15-ll-rl. This 
is a much more extreme impact of the termination shock tha n 
previously suggested in the literature (Thomp son et al.l 12001 ). 
The decelerated wind material is accumulated in a dense shell 
between the forward and reverse shocks. The pressure across 
this dense shell is nearly uniform, while the contact discon- 
tinuity between the accumulated wind matter and the dense 
layer of shock-accelerated progenitor gas is clearly visible in 
the density profiles. One should also notice that the condi- 
tions at the wind termination shock are by no means time- 
independe nt as previously ass umed in nucleosynthesis calcula- 
tions (e.g. JWanaio et al]|2002l) . Temperature and density at the 
reverse shock in model M15-ll-rl evolve, because the radial po- 
sition of the reverse shock as well as the wind properties change 
with time. The impact of the wind termination shock on the con- 
ditions in the expanding wind mass shells is better visible in 
Fig. [6] where the time-evolution of different quantities is de- 
picted as seen comoving with some selected mass shells. The 
extremely rapid decline of the temperature and density in the 
fast wind are stopped and switch over to a much slower evolu- 
tion. After the wind material has been added to the dense shell 
between the two shocks, it moves with nearly constant veloc- 
ity. Its density therefore decays approximately like p oc r 2 and 
because the gas is radiation-dominated, its temperature follows 
roughly the power law T oc r 2 ^. 

3.3. Analytic discussion of the wind termination shock 

The behavior of the wind termination shock and its effects on the 
neutrino-driven outflow can basically be understood by simple 
analytic considerations. For this purpose we consider the three 
Rankine-Hugoniot shock jump conditions for mass, momentum 
and energy flow, 



Prs Mrs — Pw^w » 

P rs +p rS M r 2 s = P w +Pw« w > 



(19) 
(20) 

(21) 



where the indices w and rs denote quantities of the wind just 
ahead of the shock and of the shocked matter just behind the 
shock, respectively. The fluid velocities m w and u rs are measured 
relative to the shock velocity, u - v - U s , P is the pressure, 
p the mass density, and cj = (s + P)/p the enthalpy per mass 
unit when s is the internal energy density of the gas. In case 
of radiation-dominated and nondegenerate conditions, one can 
write s = (s + P)/(n B k B T) for the dimensionless entropy nor- 
malized by the baryon density n B p/m B (m B is the average baryon 
mass), and therefore one gets 



s rs k B T rs - iw^B^w = ^m B (ul - Kj S ) 



(22) 



Since the wind termination shock strongly decelerates the wind, 
the postshock and preshock velocities fulfill the relation » 
u 2 s . Thus the postshock entropy is approximately given by 



k B T w 1 m B u- 
s w t-=- + - - 



k B T IS 2 k B T n 



(23) 



Again making the assumption that the gas on both sides of the 
shock is radiation dominated, the dimensionless entropy per nu- 
cleon is given by 



s = f 7 a 



(k B Tf 

«B 



(24) 



where a y — a/k^ = 2.08 x 10 49 erg~ 3 cirT 3 is related to the radia- 
tion constant a, and f y is a factor whose exact value depends on 
the temperature and thus the mixture of radiation and e + e~-pairs; 
assuming zero electron degeneracy, the corresponding range of 
values is £ < f y < y . Equation (l24l can be used to express k B T 
ahead of and behind the shock in terms of s and p. Using also 
that the densities are connected by p rs = p?p w with /3 ~ 7 for a 
strong shock and radiation-dominated conditions, one derives 



3 W 



,1/3. 



1/3 



3/4 



,4/3 
J w 

"7/3 



+ 33.5 



"w.9 

1/3 
Pw,2 



3/4 



(25) 



where a = / y a y mg/(8/3), is the wind velocity measured in 
10 9 cm s _1 , p w ,2 the wind density in 100 g cm 4 , and the numeri- 
cal value in the second expression was calculated with ft — 7 and 
f y = I. Equation (l25T l can be rewritten in terms of the wind mass 
loss rate M w and reverse shock radius R ls , using 



M w = 4nR 2 s p 



rs Pw^'w 



(26) 



and assuming that the shock velocity is negligible, and therefore 
v w = «w, which gives 



/? 2/ V /3 13/4 

■V'. , . - 1 /3 "rs "w 



jr. + 



1/3 



,4/3 
"7/3 



p 



+ 28.7 



rs,8 w,9 



M 



1/3 



3/4 



(27) 



Here /? rs g is in units of 10 cm and M w .-5 is normalized to 
10~ 5 M e . If the wind entropy is low, s w <K s, s , only the second 
terms in Eqs. d25l l and (l27T i are relevant. 

It is also possible to obtain an estimate of the reverse shock 
position from known supernova and wind parameters. In case of 
a strong shock, i.e., P, s » P w , one can derive from Eqs. $1% and 
(|20l the relation P rs m (1 - /i' l )p v/ u 2 v . Using again Eq. (|26l i for 
p w , one gets 



(28) 




P7 4;rP re 



Assuming the spherical shell between the forward shock at ra- 
dius R s » R rs and the reverse shock to have constant pressure 
and to be radiation dominated, one can make the approximation 



AnRl 



(29) 



where / exp is the fraction of the supernova explosion energy E exp 
that is present as internal energy of the gas between forward and 
reverse shock. Plugging Eq. (f29t into Eq. (f28j) one obtains 



R, 



L _ M MwUwRj 



ft } fexpE eX p 



4.14 x 10 3 



-^w^Mw.^sjo 



[km] 



(30) 
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Fig. 7. For a set of simulations with different progenitor stars, 
M10-ll-rl, M15-ll-rl, M20-ll-rl, and M25-15-r4, the different 
panels show as functions of time (from top to bottom): The radi- 
ated luminosities of v e and v e , the mean energies and rms ener- 
gies of these neutrinos (all measured at a distance of 500 km, dis- 
regarding gravitational redshift effects), the cumulative energy 
emitted in v e and v e , A£ Vc , the total energy released in neutrinos 
and antineutrinos of all flavors, AE lot , and the explosion energy 
of the models. 



The numerical value was computed by taking /3 — 7 and nor- 
malizing the shock radius to 10 10 cm, the explosion energy to 
10 51 erg, and the parameter / exp to 0. 1 . 



Equations d25l ) and ((27), evaluated with the numbers for the 
wind quantities from our numerical model, describe the entropy 
jump at the reverse shock in the simulations very well. The same 
is true for the reverse shock radius computed from Eq. d28l l with 
the pressure behind the reverse shock, P rs , taken from the hydro- 
dynamic simulations. Equation (l30l . however, does not yield a 
satisfactory agreement and for some models fails even qualita- 
tively to reproduce the behavior of the reverse shock radius as a 
function of time. The reason for this mismatch between simula- 
tions and analytic approximation is mainly the fact that the fac- 
tor /exp cannot be considered as a constant. Instead, during the 
expansion of the supernova ejecta, pdV work converts internal 
energy to kinetic energy of the matter swept up by the outgoing 
shock. Therefore, as time goes on, a smaller and smaller frac- 
tion of the explosion energy remains stored as internal energy in 
the layer between forward and reverse shock. As a consequence, 
/ exp decreases during the simulated evolution. In the first four 
seconds, / exp w 1 turns out to be a good choice, but lateron / exp 
drops monotonically to / exp » 0.25 at ten seconds. Taking this 
into account, Eq. (l30l also yields a good description of the re- 
verse shock radius as a function of time. 

3.4. Different progenitors 

The analytic discussion of the previous section, in particular 
Eqs. (l25l >. (|27"1 i. and d28T i. allow us now to understand the behav- 
ior of the wind termination shock in different progenitor stars. 
For this purpose we compare our 15 M reference model, M15- 
11-rl, with models M10-ll-rl, M20-ll-rl, and M25-15-r4, which 
are explosion simulations for 10.2, 20, and 25 M Q stars, respec- 
tively. The conditions at the inner grid boundary were chosen 
such that the models have similar explosion energies between 
roughly 1.3 B and 2 B (Table [2}. The 25 M Q star has such a big 
mass accretion rate and correspondingly high accretion luminos- 
ity that the explosion tends to become stronger than in the lower- 
mass progenitors. To lessen this effect, we reduced the boundary 
luminosities compared to the other models and chose a larger 
final radius of the inner boundary and thus of the new-born neu- 
tron star. 

The neutron star mass and radius in the 10.2, 15, and 20 M 
simulations are rather similar (Table |2]l and so are the time- 
dependent luminosities, mean energies, and energy radiated in 
v e and v e , as well as the total energy release in neutrinos of all 
flavors, AE tot , (Fig- 13- The 25 M run, however, sticks out with 
significantly higher values of all these quantities. Progenitor- 
dependent differences associated with the density structure of 
the collapsing star outside of the iron core are responsible for 
the differences in the time-dependence of the explosion energy 
seen between the 10.2, 15, and 20 M models in the lower panel 
of Fig. [7] A more massive progenitor has a higher mass accre- 
tion rate and accretion luminosity and also a larger mass in the 
gain layer. Its explosion therefore tends to be more energetic. In 
case of the models M20-ll-rl and M25-15-r4, the large binding 
energy of the outer stellar shells later on leads to a visible de- 
crease of the explosion energy from a maximum value reached 
transiently during the simulation (Fig. [7]). 

Figure [8] shows the time evolution of quantities that deter- 
mine and characterize the neutrino wind and reverse shock be- 
havior in our simulations with different progenitors. The wind 
properties (left column in Fig. [8) exhibit their well-known de- 
pendence on the neutron star mass and radius and on the neu- 
trino luminosities and mean energies. Because of the similarity 
of these quantities in case of the 10.2, 15, and 20 M models, 
only rather small differences are visible between these runs, re- 
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time [s] time [s] 

Fig. 8. Time evolution of different quantities for a set of simulations with different progenitor stars, M10-ll-rl, M15-ll-rl, M20- 
11-rl, and M25-15-r4. Shown are the baryonic mass, M\, M , and gravitational mass, M grv (Eq.|7]), neutron star radius, neutrino-wind 
expansion timescale according to Eq. ( fT5l ). wind mass-loss rate, electron fraction, and entropy per nucleon (left, from top to bottom), 
radius of the supernova shock, radius of the reverse shock, and pressure, density, temperature, and entropy per nucleon downstream 
of the reverse shock. 

vealing a slightly longer expansion timescale, lower mass-loss massive neutron star (see also Table The electron fraction 
rate, and higher entropy for model M20-ll-rl with its more shows a somewhat wider variation because of its strong sensi- 
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Fig. 9. The same as Fig. [7] but for simulations with varied inner 
boundary conditions (for clarity, however, we omitted the lines 
showing the rms energies of the radiated v e and v e and shifted the 
v e energies by 5 MeV). The models Ml 5-11 -r5, M15-ll-r6, Ml 5- 
It2-r3, and M20-13-r3 are compared with our reference 15 M 
model M15-ll-rl in order to demonstrate the influence of dif- 
ferent neutron star radii, neutron star masses, and core neutrino 
luminosities at the inner grid boundary. 



tivity to the spectral and flux differences of the v e and v e emis- 
sion. The larger neutron star mass and neutrino luminosities in 
case of the 25 M Q progenitor separate this model clearly from 
the others. They affect in particular the neutrino-wind entropy, 
which scales with the value of the neutron star mass but is only 



weakly dependent on the neutrino emission properties (Eq.[T6l>. 
Nevertheless, since the neutron star is not extremely compact 
(R ns « 1 1.5 km; Table |2]i and only moderately massive (gravita- 
tional mass M grv w 1 .66 M ), the wind entropy is never higher 
than 1 15 ks per nucleon during the 10 seconds of computed post- 
bounce evolution. The expansion timescale and mass loss rate of 
the 25 M G case are more similar to the other models because of 
a partial cancellation of their dependences on L, e, and M in 
Eqs. O and CjJ). 

The wind termination shock evolves largely differently in all 
cases (Fig. [8] right column). Obviously, the progenitor struc- 
ture has a big influence on its behavior. The supernova shock 
expands much faster in the lower-mass stars, causing a more 
rapid decline of the pressure in the shell between forward and 
reverse shock. The propagation of the forward shock and the 
time-varying conditions there are communicated inward to the 
reverse shock on the sound propagation timescale. Therefore the 
pressure just downstream of the reverse shock, f ls , as well as 
the density and temperature at this location, decrease, too. From 
Eq. d28| i it can be understood that in model M10-ll-rl the strong 
pressure reduction triggers a fast outward motion of the reverse 
shock. In the 15 M G star the increase of R rs is much less extreme, 
and in the 20 and 25 M runs the wind termination shock even 
retreats after ~2 s of initial expansion and transient stagnation. In 
these cases the decline of P re is not fast enough to compete with 
the decrease of M w and m w in the numerator of Eq. d28l ). A sim- 
ilar effect can be observed at t — Is when we change the time- 
dependence of the neutrino luminosity at the inner grid bound- 
ary. The subsequent luminosity decrease leads to the mass-loss 
rate and velocity of the wind dropping more quickly than P IS in 
all models except MlO-ll-rl, explaining why the initial expan- 
sion of the reverse shock is stopped at about this time. 

Because of the different reverse shock behavior, the den- 
sity, temperature, and entropy downstream of the reverse shock 
as functions of time show also large differences between the 
progenitors (Fig. [8). In model MlO-ll-rl the wind termination 
shock moves to radii beyond 10,000 km within little more than 
one second. During this phase the density p w behind this shock 
drops to less than 10 3 gcirT 3 and the temperature T w becomes 
lower than 10 9 K. The entropy, on the other hand, is nearly 
300 £b per nucleon after 2 seconds. In the runs for the more mas- 
sive progenitors, the density and temperature at the reverse shock 
are larger for a longer period of postbounce evolution, and the 
entropy does not reach the very high values of the 1OM simu- 
lation. The more massive the progenitor is — or, more precisely, 
the denser the shells around the iron core are — the slower prop- 
agates the shock for a given value of the explosion energy, and 
the more confined is the reverse shock. In none of the compared 
cases, however, are the conditions at the wind termination shock 
constant with time. 

3.5. Variations of inner boundary conditions 

It is clear from Eqs. $25[ , (|27| >, and (f28b that the behavior of the 
reverse shock does not only depend on the structure of the ex- 
ploding star but also on the neutrino-wind properties, in particu- 
lar the wind mass-loss rate and velocity. Since the latter increases 
with the distance from the neutron star, the radius of the reverse 
shock introduces an additional velocity dependence in Eqs. ( |2~5| >. 
63, and £Z5). 

In order to investigate the changes associated with different 
strength and time evolution of the neutrino wind, we varied the 
wind-determining parameters, i.e., the neutron star mass, radius, 
and contraction, and the core neutrino luminosities and energies 
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time [s] time [s] 

Fig. 10. The same as Fig. [8] but for simulations with varied inner boundary conditions. The models M15-ll-r5, M15-ll-r6, M15- 
It2-r3, and M20-13-r3 are compared with our reference 15 M Q model M15-ll-rl in order to demonstrate the influence of different 
neutron star radii, neutron star masses, and core neutrino luminosities at the inner grid boundary. 



as functions of time. In this section we therefore discuss the in- 
fluence of these variations of the parameters used for the inner 
boundary condition. 



The effect of the neutron star radius is visible from a com- 
parison of our reference 15 M model M15-ll-rl with models 
M15-ll-r5 and M15-ll-r6 in Fig. [10] These three simulations are 
computed with the same inner boundary condition for the neu- 
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trinos and produce neutron stars with approximately the same 
gravitational masses but final radii of 10, 13, and roughly 17 km, 
respectively (Table |2}. The neutrino luminosities radiated from 
the nascent neutron star and the energy emitted in v e and i> e as 
well as the total energy lost in neutrinos are nearly the same 
(Fig. [9]). Because of similar explosion energies, also the super- 
nova shock in the three models propagates with similar velocity 

(Fig. Hi). 

The mean neutrino energies, however, show a clear corre- 
lation with the neutron star radius: the more compact the neu- 
tron star is, the higher are the energies of the escaping v e and v e 
in Fig. [9] Also the neutrino-wind properties reveal the variation 
with the compactness of the neutron star that i s qualitatively ex- 
pected from the analytic expressions given by lOian & Wooslevl 
i 19961 see also Sect. EH Eqs. ITMISt , A larger R ns leads to a 
longer expansion timescale and thus lower wind velocity, larger 
mass-loss rate, and smaller wind entropy (see Fig. ITOb . In case 
of the mass-loss rate, however, the influence of the larger neu- 
tron star radius is partly cancelled by the lower mean neutrino 
energies (see Eq. [T8l l (and by the slightly higher neutron star 
mass of model M15-ll-r6), for which reason the differences in 
M are rather modest, in particular between models M15-ll-r5 
andM15-ll-r6. 

Qualitatively, the reverse shock exhibits the same behav- 
ior in these two models as in M15-ll-rl. While its radius R ls 
is essentially the same in models M15-ll-rl and M15-ll-r5, 
the wind termination shock, however, expands less strongly in 
model M15-ll-r6, reacting to the considerably lower wind ve- 
locity and slightly slower propagation of the supernova shock in 
this somewhat less energetic model (cf. Eqs.l28land[30l). Finally, 
the entropy of the matter decelerated in the reverse shock be- 
haves as expected from Eqs. d25l ) and d27l ) when values for the 
wind parameters and reverse shock radius are inserted into these 
equations. It is highest in model M15-ll-rl and lowest in model 
M15-ll-r6. The densities behind the reverse shock are ordered 
inversely. 

Model M15-lt2-r3 demonstrates the influence of a more 
rapid decay of the luminosities and mean energies of the radi- 
ated neutrinos after one second of postbounce evolution. This is 
associated with a reduced energy loss of the nascent neutron star 
and leads to an increase of the wind expansion timescale, a steep 
decrease of the wind mass loss rate, and a higher wind entropy 
compared to model M15-ll-rl. The reverse shock reacts to that 
by a rapid recession between about 1 s and 3 s after bounce (cf. 
Eq. l28l ) before it starts an outward motion again at later times 
when the pressure P rs drops faster than M w and m w of the wind. 
Due to the small reverse shock radius, however, u w at the shock 
is low and the entropy increase through the wind termination 
shock is modest. 

Model M20-13-r3 with a neutron star radius and neutrino 
emission properties very similar to model M15-ll-rl, but a sig- 
nificantly higher neutron star mass, reveals an even more ex- 
treme behavior. The larger neutron star mass increases the wind 
entropy, however at the same time reduces the wind mass-loss 
rate and the inverse wind expansion timescale (and thus the wind 
velocity; cf. Eqs. IT~6l - rT~8b . Moreover, the explosion energy of 
this model is very low and the supernova shock expands only 
slowly. All together forces the wind termination shock to re- 
treat as the neutrino fluxes decay, until it falls below the sonic 
point in the wind and disappears. This brings the whole region 
from the proto-neutron star surface to the outer boundary of the 
neutrino-driven outflow (which is the contact discontinuity be- 
tween shock-accelerated ejecta and neutrino-heated ejecta) in 
sonic contact, see Fig. QT| The neutrino-driven outflow is now 



only a subsonic breeze and merges with the dense shell of ejecta 
behind the outgoing supernova shock without being accelerated 
to supersonic speed. 

In order to study this phenomenon and its implications in 
more detail, we triggered the occurrence of the breeze at a much 
earlier time in model M15-ltl-r4, in which the neutrino lumi- 
nosities at the inner boundary were assumed to decay faster 
and the radiated neutrino energy is therefore even lower than 
in model M15-lt2-r3 (compare Fig. [12] with Fig. [9]). As a con- 
sequence, the neutrino wind does not have sufficient power to 
keep the wind termination shock at a large radius. In spite of 
a standard explosion energy (£ e xp ~ 1-2 bethe; Table [2] and 
Fig. [T2l and fast propagation of the supernova shock, the re- 
verse shock begins to retreat already after 1 s and disappears 
after 2.5 s (Fig. [TTJ whereas this happened only after 7 s as in 
model M15-lt2-r3 (Fig. [10). In Fig. [fe] one sees that the tran- 
sition to the subsonic breeze is accompanied by a considerable 
growth of the expansion timescales calculated from Eqs. (13[ 
and dTSt . The timescale calculated from Eq. ( TPfl i exhibits even a 
sudden increase which occurs when the wind termination shock 
has retreated so much that it is encompassed by the radial in- 
tegral of Eq. 0141 . The integral then includes shock-decelerated 
outflow which cools much more slowly. 

After about 7 seconds, however, the sound speed in the neu- 
tron star surroundings has dropped and the dense ejecta shell be- 
hind the supernova shock has moved outward sufficiently far that 
the neutrino-driven outflow can again reach supersonic veloci- 
ties, despite much less powerful acceleration than in the first two 
seconds after bounce (Fig. Iffl . This is visible also in the radial 
profiles and mass shell trajectories plotted in Fig. [13] where at 
late times (f > 8 s) the discontinuity that characterizes the pres- 
ence of a wind termination shock appears again in all quantities. 
Because of the meanwhile low wind velocity and very low mass- 
loss rate and therefore small reverse-shock radius, this shock is 
at late times much weaker than it was in the early phase. The as- 
sociated density, temperature, and entropy steps are correspond- 
ingly small (Fig.Qj}. 

During the breeze phase the outflow material is accelerated 
to a maximum velocity and then continuously decelerated again 
as it joins into the dense layer of ejecta behind the supernova 
shock. The mass-shell trajectories on the rhs. of Fi g . [T3lillustrate 
this smooth transition from the breeze expansion to the slower 
evolution when the matter is added to the dense ejecta shell. 

Models M15-ltl-r4 and M20-13-r3 demonstrate clearly that 
the wind termination shock can be a transient feature and its 
presence is very sensitive to the time-dependent conditions in 
the neutrino-driven outflow and the expansion of the dense post- 
shock shell of supernova ejecta. Simulations with a consistent 
treatment of the neutron star evolution and of the baryonic mass 
loss of the nascent neutron star are needed to make definitive pre- 
dictions of the evolution of a given progenitor star. But even then 
such predictions are handicapped by our incomplete knowledge 
of the high-density equation of state and of the corresponding 
properties and neutrino emission of forming neutron stars. 



4. Possible implications for nucleosynthesis 

Our work had the goal to investigate the fundamental aspects of 
the neutrino-driven outflow from newly born neutron stars and of 
the wind termination shock that forms by the interaction of the 
supersonic wind with the shell of slower moving ejecta behind 
the supernova shock. We intended to explore the dependence of 
this reverse shock on the wind and the explosion properties in 
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Fig. 11. Same as Fig. [3] but for model M15- 
ltl-r4. Here the explosion occurs at 0.22 s after 
bounce and because of the assumed fast sub- 
sequent decay of the neutrino luminosity, the 
reverse shock reveals a much different behav- 
ior than in case of model M15-ll-rl. It tem- 
porarily disappears when a subsonic breeze in- 



stead of a wind develops after about 2.5 s. At 
/ S 7 s the outflow expansion becomes super- 
sonic again and a wind termination shock ap- 
pears again. 




Fig. 12. Same as Fig. [4] but for model M15-ltl-r4. Compared to model M15-ll-rl, the neutrino luminosities and mean energies 
decrease faster, the total energy radiated in neutrinos and the explosion energy are lower, and the gravitational mass of the neutron 
star is larger. The breeze solutions that develop between 2.5 s and 7 s have significantly longer expansion timescales. 



different progenitor stars. The outflow conditions in our mod- 
els were found to differ significantly from the subsonic breeze 
solutions considered in many previous nucleosynthesis studies. 
Detailed network calculations are needed to analyse the conse- 
quences of these differences for the heavy-element formation in 
the outflows. Although this is beyond the scope of our present 



study, we want to add here a few speculative remarks about the 
possible implications. 

4. 1. Shocked winds vs. unshocked winds and breezes 

Two main differences of the outflows in our models are poten- 
tially relevant for the assembling of heavy nuclei, in particu- 
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Fig. 13. Left, from top to bottom: Radial profiles of the velocity, density, temperature, and entropy in model M15-ltl-r4 at different 
postbounce times. The profiles should be compared with the corresponding ones of model M15-ll-rl in Fig. [5] Right: The same 
quantities as functions of time for some mass shells ejected in the neutrino-driven outflow of model M15-ltl-r4. For comparison 
with Model M 1 5 -1 1 -rl , see Fig . 



lar the formation of r-process elements in wind phases that de- 
velop a neutron excess. Firstly, the winds accelerate to super- 
sonic speeds and therefore have shorter expansion timescales 
than breeze solutions, which by definition remain subsonic ev- 
erywhere. This affects also the temperature range from 7xl0 9 K 
down to 3xl0 9 K that is crucial for the formation of seed nu- 
clei from a particles and free nucleons through three-particle 
reactions (triple alpha and aan) and subsequent a captures. A 
faster expansion leads to less seed production and therefore an 
increase of the neutron-to-seed ratio. This increase might be sig- 
nificant for the modest values of the wind entropy obtained at the 
end of our simulations (we found up to 100-120 &b per nucleon, 
see Table [2] but 150 do not appear implausible at later times 
t > 10 s when the neutrino luminosities and mean energies have 
dropped further). Note that lTakahashi et all d 19941) . for example, 
employed breeze solutions with considerably longer expansion 
timescales in their nucleosynthesis studies. 

Secondly, while the subsonic breezes gradually slow down 
after they have reached their maximum velocity, the wind matter 
going through the reverse shock is abruptly decelerated and its 



temperature, density, and entropy are increased. In the shocked 
flow the density and temperature then continue to drop on a 
much longer timescale than during the rapid early expansion of 
the wind. If the wind deceleration happens at conditions where 
free neutrons are still present — which is the c ase at tempera- 
tures around 10 9 K in the breezes considered by Takahashi et alj 
d 1 994 ) but might be true at even lower temperatures in faster 
winds — the changed conditions may lead to a shift of the (n,y)- 
Cy,n) equilibrium and thus of the r-process path during the later 
stages of r-processing. This might alter the outcome of the nu- 
cleosynthesis compared to a freely expanding wind with its con- 
tinuously growing velocity an d compared to a gradua lly decel- 
erated breeze. The results by Takaha shi et alj (1994), see e.g. 
Fig. 2 there, suggest that the strength of the r-processing indeed 
depends on the dynamics of the outflow also at relatively low 
temperatures. In a most extreme situation (very rapid expansion, 
very high entropy) one might even imagine that neutrons remain 
unbound in the wind because their capture timescale becomes 
longer than the expansion timescale at some radius, whereas the 
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shock-decelerated environment might give all neutrons enough 
time to be absorbed by seed nuclei. 

Detailed network calculations are needed to clarify how 
strongly these differences really affect the heavy-element pro- 
duction in high-entropy supernova outflows for the conditions 
found in our models, and how the influence of these differences 
depends on the parameters that characterize the neutrino-driven 
wind and the reverse shock behavior. Such calculations being 
currently unavailable, we are not able to draw definitive conclu- 
sions here. In particular it is not clear whether indeed, and if 
so when, shock-decelerated supersonic winds can be more fa- 
vorable for a strong r-process than previously stu died breezes or 
unsho cked winds. It is interesting, however, thatjjhompson^etal] 
(2001) mentioned a strengthening of the production of third- 
peak r-process nuclei due to a slower postshock expansion and a 
shock-increased temperature at the time the r-process freeze-out 
happens. 

4.2. Progenitor trends 

We have presented simulations for a number of progenitor stars 
and have varied our prescription of the neutrino luminosities and 
of the neutron star radius and contraction with time. All of our 
choices for these conditions at the inner grid boundary were con- 
strained by our need to obtain neutrino-driven explosions (and 
reasonable values for the explosion parameters) in spherically 
symmetric models. This required the assumption of sufficiently 
large boundary luminosities or could be achieved by assuming a 
more compact neutron star. The latter case leads to an enhanced 
release of gravitational binding energy and to increased neutrino 
fluxes and thus stronger neutrino heating behind the shock. The 
considered boundary settings were also motivated by the goal 
to demonstrate the different kinds of behavior that can occur for 
the wind termination shock, depending on possible variations of 
the properties of the neutrino-driven wind, of the explosion, and 
of the progenitor conditions. Identifying clearly the influence of 
individual aspects, however, is not an easy task in simulations 
like these, where the different components that play a role are 
strongly coupled. 

Conclusions on systematic variations with the stellar progen- 
itor are not only handicapped by the limited set of investigated 
models, but in particular also by the incomplete understanding 
of the supernova explosion mechanism and of the neutron star 
equation of state. This has the consequence that neither the ex- 
plosion properties (energy, ejecta velocity, explosion timescale 
and thus mass cut) nor the properties of the compact remnant 
(e.g., its mass-radius relation and release of gravitational bind- 
ing energy) can be reliably predicted at the present time, im- 
plying that the neutrino emission as a function of time (which 
may also be strongly influenced by convection inside the nascent 
neutron star) is not well known. This in turn means that not 
only the time evolution of the neutrino-driven wind is uncer- 
tain (cf. Eqs. [T6l - [T8l ) but also the behavior and the influence of 
the wind termination shock is not definitely determined (due to 
the dependence of Eqs. [27] and [28] on the explosion and wind 
parameters). Additional complexity comes from the fact that 
multi-dimensional phenomena are important at least in some re- 
gions and during some phases of the evolution. Because of all 
these uncertainties, we concentrated here on a matter of princi- 
ple study. Future simulations that include the neutrino cooling 
of the nascent neutron star instead of making use of our inner 
boundary condition will help to link and thus to reduce the in- 
volved degrees of freedom. Nevertheless, they will still have to 
test different nuclear equations of state and will have to make 



assumptions about the onset and the energetics of the supernova 
explosion. 

In order to assess the resulting uncertainties in the prediction 
of the evolution of any particular progenitor model, we there- 
fore investigated in Sect. |3.5| models where we varied the inner 
boundary conditions within a reasonable (but maybe somewhat 
"extreme") range of possibilities. In spite of these systematic un- 
certainties, we think that our "standard" (and preferred) set of 
models in Sect. 13.41 reveals some trends that either agree with 
expectations or are likely to be confirmed by more complete and 
more detailed future studies. 

On the one hand, the neutron stars in more massive progeni- 
tors tend to be more massive, simply because of the larger stellar 
cores and the bigger mass infall rates after core bounce. This, of 
course, causes more energy release in neutrinos and higher neu- 
trino fluxes for a longer period of time. Both the larger neutron 
star mass and the differences in the neutrino emission affect the 
neutrino-driven wind. In particular they lead to a tendency to- 
wards higher wind entropies for more massive progenitor stars 
(but, astonishingly, no significant differences in the wind ex- 
pansion timescale and mass-loss rate as functions of time; see 
Fig. [S). On the other hand, more massive stars with their bigger 
and denser metal cores cause a less rapid expansion of the su- 
pernova shock during the computed phase of the explosion. The 
dense ejecta shell behind the supernova shock then forces the 
wind termination shock to follow this behavior (Fig. [8}. 

The consequence is that in less massive progenitors the re- 
verse shock reaches large radii very quickly and continues to 
travel outward. Although it can lead to a huge increase of the 
entropy in the wind, because matter is decelerated by the ter- 
mination shock from very large velocities, these effects happen 
after a short period of postbounce evolution only at a point where 
the temperature and density are already so low that the influence 
on the nucleosynthesis is probably marginal. In contrast, the re- 
verse shock in more massive stars stays at a small radius dur- 
ing the whole 10 s of computed evolution and therefore for the 
most massive progenitors the reverse shock is endangered to re- 
treat and even collapse at a later stage of the evolution when the 
neutrino-wind power has decreased. To see this happening, we 
would have to run our "standard" simulations for a significantly 
longer evolution period than just 10 seconds (with substantial 
demands on computer time because of time-step limitations and 
the need to continuously increase the numerical resolution). The 
fundamental possibility of a contraction of the reverse shock, 
however, could be demonstrated by varying the boundary con- 
ditions in some of our models such that the neutrino luminosi- 
ties and neutrino-wind power had dropped sufficiently already 
within the canonical 10 s of computed evolution (Sect. [331 ). The 
reverse shock staying at a small radius causes a more moderate 
increase of the wind entropy but still decelerates the supersonic 
wind abruptly and raises the density and temperature of the ex- 
panding matter. These discontinuous changes of the outflow con- 
ditions are time-dependent and might have interesting but so far 
unexplored consequences for the nucleosynthesis. 

We therefore conclude that the neutrino-driven outflows in 
our models are supersonic winds in which the nucleosynthesis- 
relevant conditions are sim ilar to the wind solut ions inves- 
tigated previously (e.g., by Thompso n et alj 120011) . However, 
these winds are bounded at some radius by the termination 
shock, which modifies the outflow behavior drastically. In low- 
mass supernova progenitors this reverse shock is likely to have 
nucleosynthesis-relevant consequences only during the first few 
seconds after bounce. For more massive progenitors the winds 
are abruptly decelerated, compressed, and heated by the wind 
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termination shock during a stage of their expansion (i.e. at a den- 
sity and temperature) where the r-process path might be affected. 
The combination of very rapidly expanding and ultimately su- 
personic wind and a termination shock that changes the evolu- 
tion of the outflow, has so far not been investigated in detail for 
its nucleosynthesis implications. In progenitors with very mas- 
sive and dense metal cores (roughly for stars with ^25 M ), 
which can give birth to neutron stars with a mass of ^2M , 
the wind termination shock might recede and possibly disappear 
so that the outflow develops into a subsonic breeze instead of 
being a supersonic wind. In this case the nucleosynthesis en- 
vironment becomes more similar to the condi tion s studied by 
Taka hashi etafl d 19941) . ISumivoshi etaT] d2000h or lOtsuki et all 
d2000t) where subsonic ejecta reach a maximum velocity and 
then slow down gradually to meet an imposed boundary con- 
dition at a certain radius. This boundary condition is supposed 
to account for the presence of the dense shell of ejecta that fol- 
lows the expanding supernova shock and that absorbs the later 
neutrino-driven outflow from the cooling neutron star. Our mod- 
els show that in a supernova core this "boundary" is time de- 
pendent and its consequences are not likely to be well described 
by the simple (e.g. constant) assumptions made for the ideal- 
ized neutrino-driven outflows considered in many nucleosynthe- 
sis studies. 

Our simulations therefore suggest that an analysis of the of 
the nucleosynthesis based on more detailed models of the con- 
ditions in supernova cores is very desirable. If the reverse shock 
feature has a direct relevance for the possibility of r-processing, 
its biggest impact must be expected for progenitors between 
roughly 15 M Q and roughly 25 M , where the wind termination 
shock is neither at too large radii nor disappears during the evo- 
lution phases that are most favorable for a strong r-process. 

5. Summary and conclusions 

We have presented long-time ID hydrodynamic simulations of 
neutrino-driven explosions and post-explosion outflows for pro- 
genitors with different masses. In these simulations the core 
of the shrinking nascent neutron star at neutrino optical depths 
larger than about 100 was replaced by an inner boundary whose 
contraction was prescribed and where time-dependent neutrino 
fluxes were imposed such that the neutrino-energy deposition 
around the neutron star produced explosions with a desired en- 
ergy. The time-dependence of both the boundary motion and of 
the neutrino luminosities was varied to investigate their influence 
on the high-entropy baryonic outflow and its interaction with the 
preceding supernova ejecta. Solving the Newtonian equations of 
hydrodynamics, we included approximately the effects of rel- 
ativisti c gravity by emplo ying the "effective relativistic poten- 
tial" of M arek et al.l d2006l) . This approximation yields very good 
agreement with fully relativistic calculations during the post- 
bounce accretion phase, and we found also nice consistency of 
our approach with fully relativistic solutions of stationary neu- 
trino winds. 

For the neutrino transport we performed a radial integration 
of the frequency-integrated energy and lepton number equations, 
in which the neutrino spectra were assumed to have a Fermi- 
Dirac non-equilibrium shape in the sense that the neutrino spec- 
tral temperat ure could differ fro m the temperature of the stellar 
medium (see Scheck et al. 20()(i). The neutrino luminosities and 
mean energies are thus functions of time as well as radius in our 
simulations. 

Because of the involved approximations and assumptions, 
our calculations can only be suggestive but are not suitable for 



definitive predictions of the nucleosynthesis-relevant conditions 
in dependence of the progenitor star. Our main goal was there- 
fore a matter-of -principle study of the interaction of the neutrino- 
driven baryonic outflow from the neutron star surface with the 
slower dense shell of ejecta that is accelerated by the outgoing 
supernova shock. The most important results are the following: 

- All of our models develop supersonic winds at least during 
some phases of their evolution. These winds are bounded by 
a termination shock in which the flow is abruptly deceler- 
ated and the entropy, density, and temperature of the flow 
are strongly increased. 

- The basic properties of this wind termination shock can be 
understood from simple analytic considerations using the 
shock-jump conditions at this reverse shock. The entropy of 
the shock-decelerated matter increases with the wind veloc- 
ity and is lower for high wind density. Therefore a large re- 
verse shock radius is favorable for a high entropy jump. The 
reverse shock radius increases with the mass-loss rate and 
velocity of the wind, but decreases when the pressure be- 
hind the reverse shock is high. The latter dependence links 
the behavior of the reverse shock to the propagation of the 
supernova shock and thus to the progenitor structure and the 
explosion properties. 

- The conditions at the reverse shock are progenitor-dependent 
and usually strongly time-dependent and therefore the shock 
effects are not well represented b y an outer boundary c on- 
dition with constant pressu re (e.g., Sumiyoshi et al. 200JJ) or 
constant temperature (e.g., Wanai o et al.ll2002l) . The conver- 
sion of kinetic energy to internal energy in the wind termina- 
tion shock can raise the entropy to several times the wind en- 
tropy. We found the highest values (until 10 s of postbounce 
evolution) of nearly 500 per nucleon behind the reverse 
shock (more than a factor of five increase) — but also the 
lowest temperatures (^ 10 9 K) and densities (^ 1000 g cm"" 3 ) 
— in case of the considered 10M o progenitor. In this star 
the supernova shock and the reverse shock propagate out- 
ward very rapidly. In the considered progenitors with masses 
of more than 15 M the maximum entropies are higher than 
200 ks per nucleon, corresponding to an increase of roughly 
a factor of three, with densities and temperatures behind 
the reverse shock in the first ten seconds of typically 100- 
10 4 g cirT 3 and 0.4-2 x 10 9 K, respectively. 

- When the supernova shock expands slowly (as in the case of 
very massive progenitors with big and dense metal cores) or 
the neutrino emission from the nascent neutron star decays 
rapidly and the wind power thus drops quickly, the reverse 
shock can show phases of recession and can even fall be- 
low the sonic point in the wind. The outflow then becomes 
a subsonic breeze that merges smoothly with the ejecta shell 
behind the shock without any jumps in the velocity and in the 
thermodynamic quantities. Changing conditions around the 
neutron star can lead to a re-establishment of a supersonic 
wind at later times. 

The consequences of the deceleration of the neutrino- 
driven wind or breeze by the dense shell behind the super- 
nova shock were parametrized previously in many r-process 
studies by imposing a mostly time-independ ent outer bound- 
ary condition in hydrodyna mic models (e.g., [Sumiyoshi et al.l 
2000, iTerasawa et al.l 120021) or by selecting those solutions of 
the steady-state equations that fulfill certa in conditions (e.g., 
lOtsuki et al.l2000HWanaio et al.l200U |2Q02|). Here, however, we 
have demonstrated that such prescriptions do in general not ade- 
quately account for the effects that happen when the neutrino- 
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driven outflow collides with the preceding, slower postshock 
shell. A wind termination shock can not only alter the wind en- 
tropy, density, and temperature by factors of a few, but also leads 
to a much slower expansion of the shocked outflow after its de- 
celeration. It will have to be explored in detail by nucleosynthe- 
sis calculations how the combination of the very rapid expansion 
of the supersonic winds and the changing conditions in the mat- 
ter due to the reverse shock affect heavy-element formation in 
the high-entropy supernova ejecta. 

All of our simulations exhibit the presence of proton-rich 
conditions in the neutrino-driven outflow during the first 2-3 
seconds of the explosion. Only later a neutron excess devel- 
ops in the wind matter. Our result of Y e > 0.5 in the early 
neutrino-driven wind is qualitatively in agr eement with models 
that employ Boltzmann neutrino tran sport dBuras et al.l l2006b; 
iPruet et all 120051: iFrohlich et alj |2006). We point out, however, 
that quantitatively meaningful calculations of the proton-to- 
neutron ratio require frequency-dependent neutrino transport. 
Moreover, a reliable prediction of the turning point from proton 
excess to neutron excess and of the value of Y e at late postbounce 
times is not possible without fully consistent cooling calcula- 
tions of the proto-neutron star instead of our inner grid boundary 
with prescribed time dependence of the neutrino fluxes. At best, 
our calculations can be indicative for the trends which will also 
be found in such improved ID models. 

The varying conditions in the neutrino-driven wind and the 
strong time- and progenitor-dependence of the behavior of the 
wind termination shock and of its effects on the wind raise a 
serious question: Do supernova cores provide the robust en- 
vironment for producing the extremely uniform solar-system 
like r-process abundance pattern between the Ba- and Pt-peaks 
obser ved in ultra metal-poor stars (see, e.g., ICowan & Snedenl 
2006)? In addition to the variations that are present in our mod- 
els, multi-dimensional effects lead to long-lasting a nisotropic ac- 
cretio n and at the same time directed outflows (e.g.. Sche ck et al.l 
2006). These introduce a stochastic element in the supernova 
evol ution during the first secon ds of postbounce evolution (see 
also lBurrows et al.ll2006l 120071) . The supernova ejecta in differ- 
ent directions can develop largely different conditions due to the 
strong anisotropy of the explosion mechanism and of the en- 
vironment of the forming neutron star. Downdrafts and pock- 
ets of dense, low-entropy matter in the convective shell behind 
the forward shock, for example, lag behind the overall outward 
expansion and cause a large-scale deformation of the reverse 
shock that terminates the supersonic neutrino-driven wind also 
in the multi-dimensional case. The varying radius and orienta- 
tion of the reverse shock relative to the radial direction lead to 
an angular dependence of the properties of the shocked matter. 
Therefore the amount of matter ejected with certain conditions 
(e.g., entropy, expansion timescale after passing the shock) dif- 
fers between spherically symmetric and multi-dimensional mod- 
els. A detailed hydro dynamic study in 2D is currently underway 
dArcones et al.ll2007l) . 

It is hard to see how this chaotic variability can allow for 
the robustness of environmental conditions needed for producing 
a uniform abundance pattern of high-mass r-process elements 
(even if some chunks of possible ejecta achieve to develop suit- 
able conditions of high entropy and low Y e as observed in the re- 
cent models of Bur rows et al.ll2007l) . If supernovae are the main 
sources of the high-mass nuclei beyond the A ~ 130 abundance 
peak — and a number of arguments have been made in support 
of that (e.g., ICowan & Thielemann! 12004) — a solution of this 
puzzle may be that these nuclei are produced in the later stages of 
the neutrino-driven wind, which are unaffected by the turbulent 



initial phase of the explosion. The proton-richness of the early 
supernova ejecta seen in our models in agreement with state- 
of-the-art simulations with energy-dependent neutrino transport, 
and the transition to n-richness at later times, yield support for 
this argument. The spherically symmetric supersonic winds with 
their wind termination shocks simulated in this work may be 
more characteristic of these late stages than of the early, still tur- 
bulent phases of the explosion. Previous investigations, however, 
suggest that the conditions in (unshocked) neutrino-driven out- 
flows are insufficient for a strong r-processing unless the neutron 
star i s very massive (around 2M H ) and very compact (~ 9 km; 
e.g.. lThompson et alJl200lMOtsuki et alj|2000h . 

It is currently not clear whether the supersonic outflows in 
our more "realistic" explosion models with the presence of the 
reverse shock can ease this constraint and help establishing an 
r-process favorable environment for less extreme assumptions of 
neutron star mass and radius. Maybe this is realised only in a 
subset of progenitors where the reverse shock is at a beneficial 
location. Such a requirement could single out rather discontin- 
uously progenitors with particular core and explosion proper- 
ties as favorable, while other progenitors do not develop suit- 
able conditions for a strong r-processing. More simulations, im- 
proved modeling, and detailed nucleosynthesis calculations are 
needed to explore such possibilities. 
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